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Abstract 

The approximately analytical bound state solutions of the /-wave Schrodinger equation for the 
Manning-Rosen (MR) potential are carried out by a proper approximation to the centrifugal term. 
The energy spectrum formula and normalized wave functions expressed in terms of the Jacobi 
polynomials are both obtained for the application of the Nikiforov-Uvarov (NU) method to the 
Manning-Rosen potential. To show the accuracy of our results, we calculate the eigenvalues numer- 
ically for arbitrary quantum numbers n and / with two different values of the potential parameter 
a. It is found that our results are in good agreement with the those obtained by other methods for 
short potential range, small / and a. Two special cases are investigated like the s-wave case and 
Hulthen potential case. 
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I. INTRODUCTION 



One of the important tasks of quantum mechanics is to find exact solutions of the wave 
equations (nonrelativistic and relativistic) for certain type of potentials of physical inter- 
est since they contain all the necessary information regarding the quantum system under 
consideration. For example, the exact solutions of these wave equations are only possible 
in a few simple cases such as the Coulomb, the harmonic oscillator, pseudoharmonic and 
Mie-type potentials [1-8]. For an arbitrary /-state, most quantum systems could be only 
treated by approximation methods. For the rotating Morse potential some semiclassical 
and/or numerical solutions have been obtained by using Pekeris approximation [9-13]. In 
recent years, many authors have studied the nonrelativistic and relativistic wave equations 
with certain potentials for the s- and Z-waves. The exact and approximate solutions of these 
models have been obtained analytically [10-14]. 

Many exponential-type potentials have been solved like the Morse potential [12,13,15], 
the Hulthen potential [16-19], the Poschl- Teller [20], the Woods-Saxon potential [21-23], 
the Kratzer-type potentials [12,14,24-27], the Rosen- Morse- type potentials [28,29], the 
Manning- Rosen potential [30-33], generahzed Morse potential [34] and other multiparameter 
exponential-type potentials [35] . Various methods are used to obtain the exact solutions of 
the wave equations for this type of exponential potentials. These methods include the super- 
symmetric (SUSY) and shape invariant method [19,36], the variational [37], the path integral 
approach [31], the standard methods [32,33], the asymptotic iteration method (AIM) [38], 
the exact quantization rule (EQR) [13,39,40], the hypervirial perturbation [41], the shifted 
1/N expansion (SE) [42] and the modified shifted 1/N expansion (MSE) [43], series method 
[44], smooth transformation [45], the algebraic approach [46], the perturbative treatment 
[47,48] and the Nikiforov-Uvarov (NU) method [16,17,20-26,49-51] and others. The NU 
method [51] is based on solving the second-order linear differential equation by reducing to 
a generalized equation of hypergeometric type. It has been used to solve the Schrodinger 
[14,16,20,22,48,49], Dirac [17,28,34,50], Klein-Gordon [21,24,25,50] wave equations for such 
kinds of exponential potentials. 

The NU method has shown its power in calculating the exact energy levels of all bound 
states for some solvable quantum systems. Motivated by the considerable interest in 
exponential-type potentials [12-35], we attempt to study the quantum properties of another 
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exponential-type potential proposed by Manning and Rosen (MR) [29-33] 

... . _ ( a{a - l)e-^^'' Ae'^l' \ 

where A and a are two-dimensionless parameters but the screening parameter b has di- 
mension of length and corresponds to the potential range [33]. This potential is used as a 
methematical model in the description of diatomic molecular vibrations [52,53] and it con- 
stitutes a convenient model for other physical situations. Figure 1 plots the Manning-Rosen 
potential (1) versus r for various screening distances b = 0.025, 0.050, and 0.100 considering 
the cases (a) a — 0.75 and (b) a — 1.50. It is known that for this potential the Schrodinger 
equation can be solved exactly for s-wave (i.e., I — 0) [32]. Unfortunately, for an arbitrary 
Z-states (Z 7^ 0), in which the Schrodinger equation does not admit an exact analytic solution. 
In such a case, the Schrodinger equation is solved numerically [54] or approximately using 
approximation schemes [18,50,55,56,57]. Some authors used the approximation scheme pro- 
posed by Greene and Aldrich [18] to study analytically the / ^ bound states or scattering 
states of the Schrodinger or even relativistic wave equations for MR potential [13,21]. We 
calculate and find its Z 7^ bound state energy spectrum and normalized wave functions 
[29-33]. The potential (1) may be further put in the following simple form 

V(r) = -— % , C^A,D^-A-a(a- 1), (2) 

It is also used in several branches of physics for their bound states and scattering properties. 
Its spectra have already been calculated via Schrodinger formulation [30]. In our analysis, 
we find that the potential (1) remains invariant by mapping a — > 1 — a. Further, it has a 
relative minimum value V^(ro) = ~ 4Kb2a(a-i) at tq = 61n 1 + 2"^°^^) for A/2 + a(a — l) > 
which provides 2q; > 1 + y/1 — 2A as a result of the first derivative ^|^_^^ = 0. For the case 
a — 0.75, we have the criteria imposed on the value of A is ^4 > a/2 — 3/8. For example, in 
h — jJL — l, the minimum of the potential is T^(ro) = — q;/166^(q; — 1). The second derivative 
which determines the force constants at r = Tq is given by 
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_ A^A + 2a{a - 1)] 

86%3(« _ 1)3 ■ U 



The purpose of this paper is to investigate the Z-state solution of the Schrodinger- MR prob- 
lem within the Nikiforov-Uvarov method to generate accurate energy spectrum. The solution 
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is mainly depends on replacing the orbital centrifugal term of singularity ~ [17] with 
Greene- Aldrich approximation scheme, consisting of the exponential form [16]. Figure 2 
shows the behaviour of the singular term r"^ and various approximation schemes recently 
used in Refs. [18,34,55,56]. 

The paper is organized as follows: In Section II we present the shortcuts of the NU 
method. In Section III, we derive I ^ bound state eigensolutions (energy spectrum and 
wave functions) of the MR potential by means of the NU method. In Section IV, we give 
numerical calculations for various diatomic molecules. Section V, is devoted to for two 
special cases, namely, I — and the Hulthen potential. The concluding remarks are given 
in Section VI. 



II. METHOD 



The Nikiforov-Uvarov (NU) method is based on solving the hypergeometric type second 
order differential equation [51]. Employing an appropriate coordinate transformation z — 
z{r), we may rewrite the Schrodinger equation in the following form: 

^"^^^ + ^/-^'^ + ^^-^^^ ^ °' 

where a{z) and ^(2;) are the polynomials with at most of second-degree, and t(s) is a first- 
degree polynomial. Further, using ip^i^) — (t)^{z)yn{z), Eq. (4) reduces into an equation of 
the following hypergeometric type: 

<j{z)y':,{z) + T{z)y'^{z) + Ay„(z) = 0, (5) 

where t{z) — ^^{z) -\-2'k[z) (its derivative must be negative) and A is a constant given in the 
form 

A = A„ = -nT'(z)-^^^^^<7"(^), n = 0,l,2,... (6) 

It is worthwhile to note that A or A„ are obtained from a particular solution of the form 
y{z) — yn{z) which is a polynomial of degree n. Further, yn{z) is the hypergeometric-type 
function whose polynomial solutions are given by Rodrigues relation 
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where S„ is the normahzation constant and the weight function p{z) must satisfy the con- 
dition [51] 



In order to determine the weight function given in Eq. (8), we must obtain the following 
polynomial: 

, , a'(z) —t(z) , I f a'iz) —Tiz)\^ 
<^)= ^ \ V V 2 ) -^(^) + M^)- (9) 

In principle, the expression under the square root sign in Eq. (9) can be arranged as the 
square of a polynomial. This is possible only if its discriminant is zero. In this case, an 
equation for k is obtained. After solving this equation, the obtained values of k are included 
in the NU method and here there is a relationship between A and k hy k = \ — t^'{z). After 
this point an appropriate can be calculated as the solution of the differential equation: 

- {^^^ 4>{z) = 0. (10) 

III. BOUND-STATE SOLUTIONS FOR ARBITRARY Z-STATES 

To study any quantum physical system characterized by the empirical potential given in 
Eq. (1), we solve the original SE which is given in the well known textbooks [1,2] 

2^ + ^(^)) ^(^'^' = ^^(^'^' '^)' (11) 
where the potential V{r) is taken as the MR form in (1). Using the separation method with 
the wavefunction ■0(r,^,0) = R{r)Yijn{0 , (f)) , we obtain the following radial Schrodinger 
eqauation as 



/(/+ 1) 



Rni{r) = 0, (12) 



Since the Schrodinger equation with above MR effective potential has no analjd;ical solution 

for / 7^ states, an approximation to the centrifugal term has to be made. The good 
approximation for the too singular kinetic energy term /(/ + l)r~'^ in the centrifugal barrier 
is taken as [18,33] 

- - ^, (13) 



&2 (1 — e-^/'^y 
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in a short potential range. To solve it by the present method, we need to recast Eq. (12) 
with Eq. (13) into the form of Eq. (4) by making change of the variables r — > z through 
the mapping function r = f{z) and energy transformation: 



^ = e-/^ e^^l-^J^, K.<0, (14) 



to obtain the following hypergeometric equation: 

(PR{z) ^ jl-z) dR{z) 
dz'^ z{l — z) dz 

+ 2 {-^^ + 2s^ - l{l + l)\z-[A^e^ + a{a - 1)] z^] R{z) = 0. (15) 

[z{l - z)] 

It is noted that the bound state (real) solutions of the last equation demands that 

0, when r — >■ oo, 

1, when r — > 0, 



(16) 



and thus provide the finite radial wave functions Rni{z) — )■ 0. To apply the hypergeometric 
method (NU), it is necessary to compare Eq. (15) with Eq. (4). Subsequently, the following 
value for the parameters in Eq. (4) are obtained as 

t{z) = l-z, a{z) = z-z'^, a{z) = -[A + e^ + a{a- 1)] z'^+[A + - l{l + l)\z- e^. 

(17) 

If one inserts these values of parameters into Eq. (9), with (t'{z) — 1 — 2z, the following 
linear function is achieved 

z . 1 



7r(2;) = -- ± --ya^yTa2^"+a3, (18) 

where ai = 1 + 4 [A + + a{a - 1) - k] , a2 = i {k - [A + - l{l + 1)]} and 03 = 46^. 
According to this method the expression in the square root has to be set equal to zero, that 
is, A = a\Z^ + a-iz + 03 = 0. Thus the constant k can be determined as 



A; = A-/(/ + l)±a£, a = V(l - 2a)2 + 4/(/ + 1). (19) 
In view of that, we can find four possible functions for 7r(2;) as 



z 



, z - is - %\ z, for A; = A - Z(Z + 1) + ae, 
2 I £ - (£ + I) 2; for k^ A-l{l + \)- ae. 



We must select 

k = A-l{l + l)- as, t:{z) = -- + s - (s + z, (21) 
in order to obtain the polynomial, t{z) = t{z) + 271 {z) having negative derivative as 



t{z) ^l + 2e-{2 + 2e + a)z, t'{z) = -(2 + 2e + a). (22) 
We can also write the values of A = A; + n'{z) and A„ ~ —nT'{z) — a" {z) , n = 0, 1, 2, ... 

1 



as 



A = A-Z(Z + 1) - (1 + a) 



(23) 



A„ = n(l + n + a + 2£), n = 0, 1, 2, ... (24) 

respectively. Letting A = An and solving the resulting equation for e leads to the energy 
equation 

_ (n + l)^ + Z(/ + l) + (2n + l)A-^ -1 + a 

^ " 2(n + l + A) ' ^ " ^^^^ 

from which we obtain the discrete energy spectrum formula: 



Enl = 



2^62 



(n + 1)2 + l{l + 1) + (2n + 1)A - A 



, < n,l < oo 



(26) 



2(n + l + A) 

where n denotes the radial quantum number. It is found that A remains invariant by 
mapping a —t- 1 — a, so do the bound state energies Eni- An important quantity of interest 
for the MR potential is the critical coupling constant A^, which is that value of A for which 
the binding energy of the level in question becomes zero. Furthermore, from Eq. (26), we 
have (in atomic units h — /i — Z — e — 1), 



A,^{n+1 + Ay - A(A + 1) + /(/ + 1). 



(27) 



Next, we turn to the radial wave function calculations. We use a{z) and 7i{z) in Eq (17) 
and Eq. (21) to obtain 

cf>{z) ^ z%l - z)^^\ (28) 

and weight function 

p(z) ^ z'%1 - zr^+\ (29) 



(30) 
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The functions yni{z), up to a numerical factor, are in the form of Jacobi polynomials, i.e., 
yni{z) ~ P^^^'^^'^^\l — 2z), and physically holds in the interval (0<r<oo— >-0<2;<l) 
[58]. Therefore, the radial part of the wave functions can be found by substituting Eq. (28) 
and Eq. (30) into Rni{z) = 4>{z)yni{z) as 

Rni{z) = N^iz^{l - zY^^Pt''^^'\l - 2z\ (31) 

where e and A are given in Eqs. (14) and (19) and Nni is a normahzation constant. This 
equation satisfies the requirements; Rni{z) — Oasz — 0{r^ oo) and Rni{z) = as 
z = \ (r = 0). Therefore, the wave functions, Rni{z) in Eq. (31) is valid physically in the 
closed interval z e [0, 1] or r e (0, oo). Further, the wave functions satisfy the normalization 
condition: 

oo 1 

J \Rni{r)f dr ^ 1 ^ b J z'' \Rni{z)f dz, (32) 



where Nni can be determined via 

1 

1 = bN^i J z'''-\l - z)2^+^ [Pi2^'2^+i)(l - 2z)] ' dz. (33) 



The Jacobi polynomials, Prf'''^^), can be explicitly written in two different ways [59,60]:: 

pM (^) = 2- X](-l)"-^ + + (1 - ^r-^ (1 + , (34) 

p=o VP/ P/ 

piP,^)(n _ rjn + p+l) ^fn\ r{n + p + u + r + l) /^-iV 
" ^!r(n + p + ^ + l)^ W r(r + p+l) ^ 2 ; ' ^""'^ 

where (") = ^,(^1^-,, = r(r+i)ry-r+i) • ^^^^^ ^^^^S (^4) and (35), we obtain the explicit 
expressions for Pi^^''^^~^^\l — 2z) : 

p^2s,2A+i)^^ - 2z) = (-l)"r(n + 2e+ l)T{n + 2A + 2) 



X 



-1? 



^ p\{n - py.rip + 2A + 2)r{n + 2e-p + iy 



(36) 



p(2e,2A+i)., . ^ r(rz + 2£ + 1) ^ (-l)-r(rz + 2£ + 2A + r + 2) 
^ ^ r(n + 2£ + 2A + 2) ^ r!(n-r)!r(2£ + r + l) 



(37) 
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Inserting Eqs. (36) and (37) into Eq. (33), one obtains 

1 = hNU-1) 



2, ,,„r(n + 2A + 2)r(n + 2£ + l)^ 
r(n + 2£ + 2A + 2) 

l)P+T(n + 2£ + 2A + r + 2) 



^'^'(^ ~ ^'^'^^ ~ ^^'^^^^ + + 2)r(n + 2£ - p + l)r(2£ + r + i)''*'^^' ^''^^ 



where 





Using the following integral representation of the hypergeometric function [59.60] 

1 



Re(7o) > Re(Q;o) > 0, (40) 



which gives 



1 



where 



2Fi(ao,/3o : «o + 1; l)/ao = J -^""-'(l - -2)-^«ci-2, (41) 



r(7o)r(7o - "0 - /3o) 



2Fi(ao,/3o : 7o; 1) 



r(7o - «o)r(7o - /3o) ' 



(Re(7o -ao- ^o) > 0, Re(7o) > MPo) > 0)- (42) 

For the present case, with the aid of Eq. (40), when aQ = n + 2e + r— p, /3q — —p — 2A — 2, 
and 7o = CKo + 1 are substituted into Eq. (41), we obtain 

^ 2^1 jao, : 7o; 1) ^ r(n + 2£ + r - p + l)r(p + 2A + 3) 
"'^^'"■^ ao (n + 2£ + r-p)r(n + 2£ + r + 2A + 3)' ^ 



Finally, we obtain 



1 _ ^j^2 ( ^.. r(» + 2A + 2)r(n + 2c + l)^ 
^"^^"'^ r(n + 2£ + 2A + 2) 



(-l)P+T(n + 2£ + r - p + l)(p + 2A + 2) 

^ p!r!(n-p)!(n-r)!r(n + 2£-p+l)r(2£ + r + l)(n + 2£ + r + 2A + 2)' ^ 
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which gives 



Nni = (45) 

where 

sin) - hi iNn r(n + 2A + 2)r(n + 2. + ir 

^ (-irT(n + 2£ + r - p + l)(p + 2A + 2) 

^ p!r!(n-p)!(n-r)!r(n + 2£-p+l)r(2£ + r + l)(n + 2£ + r + 2A + 2)' ^ 



IV. NUMERICAL RESULTS 



To show the accuracy of our results, we calculate the energy eigenvalues for various n and 
/ quantum numbers with two different values of the parameters a. Its shown in Table 1, the 
present approximately numerical results are not in a good agreement when long potential 
range (small values of parameter h). The energy eigenvalues for short potential range (large 
values of parameter h) are in agreement with the other authors. The energy spectra for 
various diatomic molecules like HCl, CH, LiH and CO are presented in Tables 2 and 3. 
These results are relevant to atomic physics [61-64], molecular physics [65,66] and chemical 
physics [67,68], etc. 



V. DISCUSSIONS 



In this work, we have utilized the hypergeometric method and solved the radial SE for 
the M-R model potential with the angular momentum I ^ states. We have derived the 
binding energy spectra in Eq. (26) and their corresponding wave functions in Eq. (31). 

Let us study special cases. We have shown that for a = (1), the present solution reduces 
to the one of the Hulthen potential [16,19,57]: 

yW(^) ^ -Hj^^, ^0 = Ze^b, 5 = b'' (47) 

where Ze^ is the potential strength parameter and 5 is the screening parameter and b is 
the range of potential. We note also that it is possible to recover the Yukawa potential by 
letting 6 — )> oo and Vq = Ze^/b. If the potential is used for atoms, the Z is identified with 
the atomic number. This can be achieved by setting A = /, hence, the energy for / ^ 



10 



states 



[A-{n + l + lffn^ 
8/x62(n + / + 1)^ 



nl 



, < n, Z < oo. 



and for s-wave (/ = 0) states 

En 



[>l-(n + l)f fi^ 
^ ^2~' 0<n<oo 



(48) 



(49) 



8/i62(n + l)^ 

Essentially, these results coincide with those obtained by the Feynman integral method 

[31,56] and the standard way [32,33], respectively. Furthermore, if taking h = 1/5 and 
identifying as Ze^5, we are able to obtain 



2n2 



2h^ 



-I 2 



which coincides with those of Refs. [16,19]. Further, we have (in atomic units % = ii = Z 

e=l) 

1 (n + / + l) 



(50) 



E. 



nl 



n + l + l 



(51) 



which coincides with Refs. [16,33]. 

The corresponding radial wave functions are expressed as 



where 



1^ 



in + l + 1] 



< n.l < oo. 



(52) 



(53) 



n + l + 1 2Ze2/x 

which coincides for the ground state with that given in Eq. (6) by Goniil et al. [18]. In 
addition, for 5r <^ 1 (i.e., r/b <^ 1), the Hulthen potential turns to become a Coulomb 
potential: V{r) — —Ze^/r with energy levels and wave functions: 



E, 



nl 



{n + l + 1)^ 



, n = 0,l,2,.. 



2*2 



z'h 



■So 



2/iao 



2 ' 



(54) 



where Eq — 13.6 eV and oq is Bohr radius for the Hydrogen atom. The wave functions are 

L h [n + l + l)_ 
which coincide with Refs. [3,16,22]. 
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VI. CONCLUSIONS AND OUTLOOK 



In this work approximately analytical bound states for the Z-wave Schrodinger equation- 
with the MR potential have been presented by making a proper approximation to the too 
singular orbital centrifugal term ~ r~^. The normalized radial wave functions of Z-wave 
bound states associated with the MR potential are obtained. The approach enables one 
to find the /-dependent solutions and the corresponding energy eigenvalues for different 
screening parameters of the MR potential. 

We have shown that for a = 0, 1, the present solution reduces to the one of the Hulthen 
potential. We note that it is possible to recover the Yukawa potential by letting 6 — > oo and 
Vq — Ze'^/b. The Hulthen potential behaves like the Coulomb potential near the origin (i.e., 
r — >■ 0) Vc{r) = —Ze'^/r but decreases exponentially in the asymptotic region when r S> 0, 
so its capacity for bound states is smaller than the Coulomb potential [16]. Obviously, the 
results are in good agreement with those obtained by other methods for short potential 
range, small a and 1. We have also studied two special cases for / = 0, / 7^ and Hulthen 
potential. The results we have ended up show that the NU method constitute a rehable 
alternative way in solving the exponential potentials. We have also found that the criteria 
for the choice of parameter A requires that A satisfies the inequality Vl — < 2a — 1. 
This means that for real bound state solutions A should be chosen properly in our numerical 
calculations. 

A slight difference in the approximations of the numerical energy spectrum of Schrodinger- 
MR problem is found in Refs. [55,56] and present work since the approximation schemes 
are different by a small shift In our recent work [17], we have found that the physical 

quantities like the energy spectrum are critically dependent on the behavior of the system 
near the singularity (r = 0). That is why, for example, the energy spectrum depends 
strongly on the angular momentum /, which results from the singularity of the orbital 
term, even for high excited states. It is found that the orbital term is too singular, then 
the validity of all such approximations is limited only to very few of the lowest energy states. 
In this case, to extend accuracy to higher energy states one may attempt to utilize the full 
advantage of the unique features of Schrodinger equation. Therefore, it is more fruitful to 
perform the analytic approximation of the less singularity rather than the too singular 
term which makes it possible to extend the validity of the results to higher excitation 
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levels giving better analytic approximation for a wider energy spectrum [69] . 
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FIG. 1: Variation of MR potential as function of separation distance r taking various values for 
the screening parameter b when (a) a = 0.75 and (b) a = 1.50. 



FIG. 2: A plot of the variation of the singular orbital term (dotted-solid line) with the 

approximations of (a) Ref. 34 (dash line), the conventional Greene- Aldrich of Ref. 18 (dash-dot 
line) and improved [55,56] (solid line) replacing the term 1/r^ with respect to r where 6 = 0.1 
fm~^, and (b) the improved approximation [55] with various shifting constants. 
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TABLE I: Energies (in atomic units) of different n and / states and for a = 0.75 and a = 
A = 2b. 

a = 0.75 a = 1.5 

states 1/6 Present QD [33] LSI [54] Present QD [33] LS [54] 
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TABLE II: Energy spectrum of HCl and CH (in eV) for different states where Tic = 1973.29 eV 
A°, iijjci = 0.9801045 amu, = 0.929931 amu and A = 2b. 
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TABLE III: Energy spectrum of LiH and CO (in eV) for different states where He = 1973.29 eV 
A°, ^j^ijj = 0.8801221 amu, fico = 6.8606719 amu and A = 2b. 
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